CERN-TH/96-301 



hep-ph/96 11227 



Neutral Heavy Leptons and Electroweak Baryogenesis 

P. Hernandez and N. Rius Q 

Theory Division, CERN 
CH-1211 Geneva 23, Switzerland 



Abstract 

We investigate the possibility that baryogenesis occurs during the weak phase transition in a 
minimal extension of the Standard Model which contains extra neutral leptons and conserves total 
lepton number. The necessary CP-violating phases appear in the leptonic Yukawa couplings. We 
compute the CP-asymmetries in both the neutral and the charged lepton fluxes reflected in the 
unbroken phase. Using present experimental bounds on the mixing angles and Standard Model 
estimates for the parameters related to the scalar potential, we conclude that it seems unlikely to 
produce the observed baryon to entropy ratio within this kind of models. However, we comment 
on the possibility that the constraints on the mixings might be naturally relaxed due to small 
finite temperature effects. 
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1 Introduction 



The baryon number to entropy ratio in the observed part of the Universe is required to be Ub/s ~ 
(4-6) x 10~ n by nucleosynthesis constraints |]T]. In 1967, Sakharov [[| established the three basic 
requirements for obtaining this baryon asymmetry as a result of particle interactions in the early 
universe: a) Baryon number violation, b) C and CP violation, c) departure from thermal equilibrium. 
These conditions may be fulfilled at weak scale temperatures ||, if the electroweak phase transition is 
first order. 

In a strongly first order electroweak transition, bubbles of the true ground state (broken phase) 
nucleate and expand until they fill the Universe; local departure from thermal equilibrium occurs in 
the vicinity of the expanding bubble walls. The other two Sakharov conditions are also satisfied, since 
C and CP are known to be violated by the electroweak interactions and anomalous baryon number 
violation is fast at high temperatures in the symmetric phase. As a bubble expands, particles in 
the unbroken phase will reflect off the advancing wall. CP-violating interactions result in a different 
reflection probability for fermions with a given chirality and the corresponding antifermions, leading 
to a CP asymmetry in the reflected chiral number flux |4|. In the symmetric phase, anomalous 
B + L violating interactions are in thermal equilibrium and the reflected current induces a net baryon 
number. An important survival requirement for the produced baryon asymmetry is that the sphaleron 
processes inside the bubble are slow enough and this in turn is directly related to the strength of the 
phase transition. 

In principle the Standard Model (SM) contains all the necessary ingredients for electroweak baryo- 
genesis, but it has two problems: the CP asymmetry induced by the Kobayashi-Maskawa phase is far 
too small to account for the observed n B /s ratio [§, ||, and the phase transition appears too weakly 
first order for the Higgs mass experimentally allowed However, these two problems may be absent 
in several simple extensions of the SM, which contain additional sources of CP violation and more 
scalars than the SM. The larger parameter space in the scalar sector allows for a stronger first-order 
phase transition without such a light Higgs ||. Several such possibilities have been considered in the 
literature: two Higgs models with a strong CP phase |P|— 1|12||, heavy Majorana neutrinos [13|, and 
supersymmetric models WA . 



In the present paper, we consider models with an extended lepton sector, which conserves total 
lepton number. The model provides a viable alternative to the see-saw mechanism for explaining the 
lightness of the known neutrinos, in those extensions of the SM where there are no scalars carrying 
lepton number, which could generate Majorana masses. They arise in several contexts, such as GUTs 



15| and E(6) superstring-inspired models |l6fl. Similar patterns of lepton masses have also been 



obtained in the context of models of Extended Technicolor with a GIM mechanism [17| . 

The relevant features for baryogenesis are twofold. First, the lepton Yukawa interactions contain 
additional CP-violating phases, which can lead to a much larger CP asymmetry than the CKM phase 
in the SM. In contrast to the models with Majorana neutrinos considered in |4j], the CP- violating effects 
in this case are not suppressed by the light neutrino masses. Second, the presence of an additional 
singlet scalar may help in getting a stronger first-order phase transition. 

An interesting issue that may be relevant in this type of models is whether finite temperature 
corrections can produce an enhancement of the CP asymmetry, as was found in the first detailed 
calculation of this quantity in the SM |18| []. In |5j || it was shown that this enhancement disappears 

lr ro be more precise, the typical suppression with quark masses expected in a flavour-blind CP-violating process was 
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in the SM when one properly includes the incoherence effects induced by the interaction of the quarks 
with the plasma. However, the question remained that if the particles involved were much more weakly 
interacting, as leptons instead of quarks, maybe this enhancement would be at work. 

In section 2 we describe the model, and the order of magnitude estimates of the CP-violating 
asymmetries are obtained in section 3. We find that the leading effects are different for the reflection 
of the neutral and the charged leptons. Naively the later is smaller since it is suppressed by the charged 
lepton masses; however, an enhancement of the type of ref. fl8| could imply no such suppression. In 
section 4 we compute the contribution to the asymmetry due to the reflection of the neutral leptons, 
which turns out to be the leading effect, as expected. In close analogy with the SM case, we consider 
in section 5 the lepton asymmetry generated by the charged lepton reflection on the bubble wall. As 
we will see, no enhancement with respect to the naive estimate is found. In section 6 we compute the 
baryon number induced by the CP asymmetries in the neutral sector, and we conclude in section 7. 



2 The Model 



The phenomenology of this type of models has been extensively studied in [19|[17|. Here we briefly 
describe the essential features relevant for baryogenesis. 

The gauge group is the standard SU(2) x U(l), with minimal quark sector. The lepton sector is 
extended with two electroweak singlet two-component leptons in each generation, i.e., 



n = [ e f ) . u R , si. (1) 

Unlike the minimal standard model, total lepton number conservation is not an automatic symmetry. 
It has to be imposed, and it restricts the form of the Yukawa terms that lead to the neutral fermion 
masses, while the Yukawa terms involving charged leptons are completely standard: 

C Y = ^ L fHe R + ^ L f D Hv R + s L f s au R + h.c. (2) 

where /j are the Yukawa matrices, H is the standard Higgs doublet and a is a new singlet scalar field. 
Due to the presence of a, the weak phase transition can be quite strongly first order for a significant 
range of parameters |§. 

For simplicity, we will assume that the singlet a acquires a vacuum expectation value at the same 
scale as the standard Higgs doublet, and the theory undergoes a single phase transition at a critical 
temperature near the weak scale. Although this does not need to be the case, we expect that the CP 
asymmetry will not depend much on this choice (at least if power counting arguments give a correct 
estimate). Then, in the broken phase the lepton mass terms are 

e R me L + vrDv l + u R Ss L + h.c, (3) 

where 

m = fv D = f D v S = flu (4) 
v and u being the vevs of the doublet and singlet scalar fields, respectively, 
not found, while the suppression with the CP-conserving angles was explicit. 
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The structure of the neutral sector mass matrix (|3|) ensures the existence of three massless Weyl 
neutrinos z/°, regardless of the relative value of D and S. The other six Weyl fermions pair up into 
three heavy neutral Dirac fermions n, whose masses are essentially determined by S. Note that in this 
model the ratio D/S is expected to be less constrained than the corresponding parameter of any model 
that invokes the see-saw mechanism to understand the smallness of neutrino masses. This is because 
in such models the neutrino masses are only suppressed by D/S, which is therefore very constrained. 
In the present case, this ratio is not related to the light neutrino masses. Nevertheless, we will see in 
the next section that it can also be constrained from the strong bounds on charged lepton mixing. 

The mass matrix can be diagonalized by multiplying on the right and the left by unitary matrices: 

(5) 
with 

where u° are the massless neutrinos and n are the neutral heavy leptons (NHL). 

The unitary matrix V diagonalizes DD* + SS^ to give M 2 . The unitary matrix U can be written 

as 

U=( Kl Kh ) (7) 
^ K SL K SH ) [i) 

where K L , K SH ~ 1 - 0[(D/S) 2 } and K H , K SL ~ 0(D/S). 

Using (0) we get the following form for the charged current leptonic weak interaction: 

£ cc = 9 - W» £ E W^L \K L u° + K H n L ] + h.c. (8) 

One can see that the charged current coupling of the mass eigenstates charged leptons to the massless 
as well as the heavy neutrinos is non-trivial, making possible the violation of individual lepton numbers 
L e , L ^ and L T . 

Due to the admixture of fermions of different weak isospin, there is no GIM mechanism in the 
neutral fermion couplings to the Z boson, which are given by 

C nc = -£-Z» N altl LP aP N p , (9) 

a,/3 

where N a = (i/®,n a ), and P = K^Kl + K\jKh is in general a non-diagonal projection matrix. The 
neutral couplings involving the massless neutrinos are diagonal but flavour-dependent. 

It has been shown [^(| that for n generations the total number of physical parameters describing 
the Yukawa sector is n 2 angles and (n — l) 2 phases. Thus, for three families there are four independent 
CP- violating phases. If the charged lepton Yukawas are neglected, it is easy to show, using the method 



of ref. [21], that the number of phases is (n — l)(n — 2)/2. So one CP-violating phase still remains in 
this limit and consequently the associated CP invariant is not suppressed by the charged lepton masses. 
In contrast, the CP invariants corresponding to the other three phases are necessarily suppressed by 
the small differences of charged lepton masses. 
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3 Order of Magnitude Estimates 



The observable CP asymmetry results from the interference of pure CP-violating phases with CP- 
even phases, equal for particles and antiparticles. These are the reflection coefficients, which become 
complex when the particle energy is smaller than its mass in the true vacuum. Unremovable CP-odd 
phases appear in the mass matrices due to either 

a) CP- violating interactions in the thermal loops that correct the dispersion relations of the particles 
propagating in the plasma |D|] ; 

b) non-trivial space-time dependence of the scalar vevs inside the bubble wall (for more than one 
Higgs field), which induces space dependent CP- violating phases. These phases cannot be rotated away 
at two adjacent points, x and x + dx, by the same set of unitary transformations, i.e. U~ l U x+ dx 7^ 1 
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Whenever mechanism b) is present, it generically will dominate over a), since in a) there are 
suppression factors coming from loops (l/Air). Mechanism b) is the one that generates the baryon 
asymmetry in all the extensions of the SM proposed in the literature for electroweak baryogenesis. In 
contrast, in the SM the quark mass matrix has only an overall dependence on the Higgs vev and can be 
diagonalized by space-independent unitary matrices; hence the CP asymmetry can only be generated 
through mechanism a). 

In the model considered here we have to distinguish between the charged and the neutral sectors. 
Charged leptons get their masses only from the doublet scalar vev, so the situation is completely 
analogous to the SM: CP-violating phases appear in the thermal corrections to the dispersion relations. 
In the neutral sector the situation is different because the mass matrix has a non-trivial dependence 
on both singlet and doublet scalar vevs. Since generically this ratio is not constant within the wall, 
mechanism b) is also present. 

The size of the leading CP asymmetries in the reflection of both charged and neutral leptons can be 
estimated by simple power counting arguments. To do so, we construct a measure of the CP violation, 
invariant under flavour and phase redefinitions of the lepton fields, i.e. under transformations of the 
type 

e R -> Ve R 

vr -> Wv R (10) 
s L -> Xs L . 

One can show that the following expression is invariant under such transformations, and vanishes if 
CP is conserved |2Cj: 



ImTr 



D^Dm)mD ] SS ] D 



[ID 



Notice that this effect cannot be tree level in the reflection amplitude, since it involves the couplings 
of both charged and neutral leptons. Therefore it is typically down by loop factors (l/47r). 
This invariant is given by 

Sep = E M 2 a M 2 b {Ml - M 2 b ) £K 2 - m))lm{K Hia K* Hja K* Hih K Hjb ), (12) 

a<b i<j 

where Mi and m a are the NHL and charged lepton masses, respectively. The natural scale in the 
problem is of the order of the electroweak phase transition temperature, T ~ 100 GeV. Therefore, 
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to obtain a dimensionless quantity 5q P should be divided by T 8 . Typically Mj ~ T, but the small 
charged lepton masses give a suppression of order (m T /T) 2 ~ 10~ 4 at least. We expect from eq. (|i~2D 
that the leading CP asymmetry in the reflection of charged leptons will appear at fourth order in the 
mixing K H ~ 0(D/S). 

On the other hand, the leading effect coming from the neutral sector appears at tree level. The 
leading CP measure, invariant under the transformations of ([II]) involving only the neutral fields, is 
given by 



ImTr 



DD*S&(DDi) 2 (SSty 



which gives 



S 3 cp 



M 2 M 2 M 2 (M 2 - M|)(M 1 2 - M 3 2 )(M 2 2 - M 3 2 )Im \{K ] H K H ) l2 {K ] H K H ) 2 ^K ] H K H ) u 



(13) 



(14) 



In this case, the asymmetry appears at sixth order in the mixing, 0[(D / S) e }. To obtain a dimensionless 
quantity, we should divide S^ P by T 12 but there is no suppression in the masses here since Mj ~ T. 
There is also an additional contribution of the form ([12]) due to loop corrections, but it would be 
suppressed at least by (m T /T) 2 ~ 10~ 4 and by loop factors, which considering the experimental 
bounds, is a larger suppression than the extra (D/S) 2 . 

If we assume that the asymptotic value of the ratio of scalar vevs is the same as at zero temperature 
v(T)/u(T) = v (0)/w(0), there are quite strong experimental bounds on the elements of the submatrix 
K ] H K H . These bounds depend on the NHL mass @§. For 3 GeV < M < M z the strongest 
limits come from LEP and they are very stringent Q: \hn(K H iaK* H j a K* Hib K H jb)\ < 3 x (10~ 9 -10~ 7 ). If 
M > Mz, there are low-energy constraints that arise both from the non-observation of lepton flavour 
violation and from universality, as well as limits from the invisible width of the Z boson ||23|| . The 
limits are slightly weaker for the mixing of the third family with any of the first two, 



\lm(K Hia K* H3a K* mb K H3b ) 



< 5 x 1(T 5 (10" 



(15) 



The first number corresponds to the so-called 'joint' bounds in ref. [2~3fl, for which cancellations among 
the different possible fermion mixings are allowed, while the number in brackets corresponds to the 
'single' limits, obtained when the remaining mixing parameters are set to zero. 



For the invariant ([13]), we get 



Im 



{K\jK h ) 12 {K^ h K h ) 2 ^ h K h )J < 1(T 7 (3 x 10 



(16) 



Based on the bounds (|H|) and taking into account the loop factor expected in that case, plus the further 
suppression in the charged lepton masses, typically of order (m T /M) 2 ~ 10~ 4 , the CP asymmetry in 
the reflection of charged leptons (|l"2"D is expected to be too small to generate a significant baryon 
asymmetry. However, a similar enhancement as the one found in |TB| could imply that there is no 
suppression coming from the lepton mass, in which case the effect could be important. This is the 
reason why we decided to do a detailed calculation in this case. 

In the case of the reflection of the neutral leptons, the bound (|T6|), together with the fact that there 
is no power suppression in the light masses or loop factors, implies that the effect could be of roughly 
the right order of magnitude. 



2 These bounds have been obtained using inequalities of the form llm^HiaK^^K^^Knjb)] < 
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4 CP Asymmetry in the Neutral Sector 



In this section we will compute the CP asymmetry in the number current of Ul that get reflected into 
the unbroken phase. This asymmetry in the rest frame of the wall is simply given by 



• _ -trans , -trans , -ref 
JCP — Js b L ->v£ ' Jv\-^v1 "T Jv^vl 



J 0^ EUW-^iV* M + M ')i 2 - \ A k^i^ M + m!)\ 2 } v ^f\{v) 

+ / 7^ E[I^ M (-P.ft)| 2 - \AU^t(-p z ,P z )\ 2 ] |/o (P), (17) 
with being the sum over flavours and 

^o(p) = e (\p\+v w Pz )/t + p ( 18 ) 



^ e (b1-«)/T + 1 > #(p0 e(v ^j5^-« IB p i )/T + 1 ' (19) 

being the thermal distributions of the different particles in the unbroken (u) and broken (b) phases as 
seen from the rest frame of the wall; v w is the wall velocity, which is estimated to be v w ~ 0.1-0.4 in 



the SM |2§ 



In the present case, thermal corrections to the propagation are negligible. The thermal masses are 
< 0.25Mj for the heavy leptons and thus considerably smaller than the energies at which the effect will 
be significant, u > min{Mi}. Furthermore, the mean free path of these weakly interacting particles is 
expected to be large compared both to the expected width of the bubble wall and to the reflection time 
of the leptons ~ Mf ; in the scattering with the wall the neutral leptons will therefore be assumed to 
be free. The transmission and reflection amplitudes will thus be computed at zero temperature, using 
LSZ reduction formulae in terms of the propagator in the presence of the wall: 

A = J d 4 x J d 4 y e - iq * x e iq ? y u(q f )(id - m)S(y,x)(-id - m)u(ft) 

= (2vr) 3 %* - qf ) S(qj - qf) 5(E f - E t ) A{q*, q}), (20) 

with 

S(y,x) = (0\T[*(y)*(x)]\0). (21) 

An analogous expression holds for antiparticles. The spinors in formula fl2"UD are on-shell and normalized 
to unit flux in the z direction, i.e. 

M7 2 M = 1. (22) 
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Since the potential created by the bubble wall is only dependent on the coordinate z, momenta in 
the x and y directions are conserved. The transmission and reflection amplitudes only depend on the 
momenta in the z direction and can be computed in a much simpler way by first boosting to a frame 
where q x ,q y = 0. With the proper normalization chosen for the spinors (|22|) , the amplitude in the 
boosted frame is simply given by (^), with the propagator and incoming and outcoming momenta 
substituted by the boosted ones. 

We can further simplify the expression for jcp by using CPT symmetry and unitarity constraints, 
which imply 

E \ A lil-> u il\ 2 = 1 ~ E l^^^l - E ( 23 ) 

i " " i i 

and substituting (|23|) in eq. (|17|): 



+ / 7^ E(K^«| 2 - |(/f b') - /oW)- (24) 

Finally, by expanding the Fermi distributions for small wall velocities: 



t y (2tt) 



/p 2 — M 2 — v? \ r> 2 

+E r - 1 )(I4^I 2 - l^k^l 2 ) %Mp)[i - Mp)1, (25) 



where = l/(e'^ //T + 1) is the unboosted Fermi distribution. 

In order to compute the amplitudes in eq. (p0| ) we would need the exact propagator in the presence 
of the wall. The potential created by the wall in the weak basis (z/ L , sl, u r ) is 



M(z) 



I V{z)K H M 

U{z)K SH M | , (26) 

V V{z)MK\j U(z)MK ] hs 



where V(z) = v(z)/v is the ratio of the vevs of the doublet scalar H in the wall and the asymptotic 
vev in the broken phase; U(z) is the ratio corresponding to the singlet a field and M = {Mi) is the 
diagonal mass matrix of the Dirac neutrinos. 



The simplest approach would be to do perturbation theory in M(z) || |14 |, which is effectively 
an expansion in M{z)/uj. Although this approximation makes the calculation much simpler, it is not 
justified since the region of interest is always to ~ M*. Instead, we will perturb in the mixing, that is 
in Kh = 0(D/S) and Khs — 1 = 0[{D / 'S) 2 }. We can write the mass matrix as 

M(z) = M (z) +5M(z), (27) 



with 





M (z) = | U(z)M | , SM(z) = M{z) - M (z). (28) 
U{z)M 
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Our strategy will be to solve the scattering problem with potential Mq(z) exactly and perturb 
only in 8M{z). This can be done for several forms of U(z). In order to simplify the problem we 
consider here the thin wall approximation for the singlet field i.e., U(z) = 6(z). The result will be 
more important for a singlet width as different from the doublet width as possible, so we will keep the 
singlet width to its minimum value and vary the doublet one. There is no reason to expect that any 
other choice would give very different results. The perturbation in 5M then gives 

S(x 2 , Xl ) = [Hdzi S®(x 2 , z l )SM(z 1 )S^(z 1 ,z 2 )5M{z 2 ) . . . S^^x,), (29) 

where the integration is done over all 00,00) and is the exact propagator in the potential 
Mq(z). It is a matrix with spin structure 



S (0) 



( \ 



sffi sffi 

(0) Q (0) 



(30) 



\ Srl Srr J 



We can then use the zero temperature propagator in the presence of a thin wall that has been computed 
in [P5[| . The formulae are given in appendix A. 

We approximate the doublet field wall profile, V(z), as 

z < 

V(z) = { 5 H l z 0<z<5 H (31) 

1 z > 5 H 

so that the wall thickness is parametrized by 8jj. We expect that this simple form is enough to give a 
reasonable estimate of the CP asymmetry. 

The calculation is straightforward. In the case of three families, the result turns out to be non- 
zero at sixth order in (D/S) as expected from the invariant ([□!)■ The contribution coming from the 
reflection is 

Y,(\ A > R ^t I 2 - antiparticles) = J ijk F r {M t , Mj, M k ), (32) 

and the transmission one 

E<K^-^ I' ~ antiparticles) = J ijk F\M U M jt M fc ), (33) 



'•j 



where 



J ljk = \m[{K\jK H)ij{K\jK}j) ki {K\jK h 



F r {M iy Mj, M k ) = 2M}MjMZ{Im[A' k A ij } - 2Im[7*^-)]} 

^M„M„M k) = -MfMfM! ^ |; + -;J l m {gg- + 2/;;^-) I , (34) 
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Figure 1: eiji ! F t (M i , Mj, M k ) as a function of p z for NHL masses Mj = (.6, .8, 2.) and S H = .1 (solid), 
dtf 1 = .2 (dashed-dotted) and Sjj 1 = .3 (dashed). All in units of the temperature. 



The integrals Ix a , . . . , I 5a are defined in appendix B. 

It can be easily checked that whenever two masses are degenerate the result vanishes. We have also 
checked that in the thin wall approximation, i.e. 8h —>■ 0, the effect disappears, as it should happen 
since in this limit both the singlet and doublet wall profiles become the same. 

As we saw in the previous section, for three generations there is only one CP-violating phase and 
Jijk — J123 tijk- In this case, the phase Ji 23 factorizes in the asymmetries of eqs. ( |3~2"1) and d33f ). We 
find that the contribution from the transmission amplitude is the dominant one, while the reflection 
amounts to a small correction. In figs. 1 and 2 we plot €ijkF t (M i , Mj, M^) as a function of p z , for 
different values of the NHL masses and Sh- We will consider masses of order ~ T, because for heavier 
neutral leptons the current will be strongly suppressed by the Fermi distribution in eq. ([251). 

Figure 1 shows the dependence on Sh for a fixed NHL spectrum. We find that the effect is more 
important when the mass differences of at least two NHLs are (Mj — Mj) < 5 H . When all mass 
differences are larger than 8 H , the asymmetry oscillates rapidly (we expect the oscillation period to 
be related to 5h) and the integrated result is suppressed. Also we observe that the effect is smaller for 
smaller 5h, in agreement with the fact that it vanishes in the limit 5h —> 0. Thus we expect that the 
largest effect will occur for Sh, satisfying (Mj — Mj) < S H , as large as possible. 

In fig. 2 we fix 5h and keep M 2 — M 1 < Sjj 1 , while varying the masses M 1 and M 3 . The peak 
appears after the second threshold (i.e. p z > M 2 ); thus, as M 2 changes, the position of the peak moves 
accordingly. Nevertheless, the integrated result is not very sensitive to the particular values of the 
NHL masses, provided the relation M 2 — Mi < 5jj l is satisfied. 

The current jcp we obtain for generic values of the masses, with the only restriction that Mj ~ T 



and 




(35) 
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Figure 2: €ijkF t (M i , Mj, M&) as a function of p z for 5^ = .3 and NHL masses: i) (.8, 1.1, 1.4) solid 
line, ii) (.8, 1.1, 2.) dashed-dotted, iii) (.6, .8, 1.) dashed and iv) (.6, .8, 2.) dotted. All in units of 
the temperature. 

and Mi — Mj ~ djj 1 for at least two NHL, is typically (0.5-1) x 10~ 2 Ji 2 3 v w . 

5 CP Asymmetry in the Charged Sector 

In this section we compute the CP asymmetry in the flux of charged left-handed leptons, reflected 
in the unbroken phase. The charged lepton mass matrix has just an overall dependence on the doublet 
scalar vev (as occurs in the SM); therefore only mechanism a) as defined in section 3 is present in this 
case. 

The calculation of the charged lepton CP asymmetry is completely analogous to the computation 
done for quarks in the SM |18|, |5|, |(| , and we refer the reader to these works for further details. Contrary 
to the case of the NHL, the one-loop thermal corrections are much larger than the tree-level masses of 
the charged leptons. The resummation of the thermal self-energies considerably modifies the dispersion 
relations and the correct asymptotic states are now quasi-particles. 

Following the notation of ref. 



26], the thermal one-loop contribution to the charged lepton self- 



energy in the broken phase can be written as 



Re(E(fc)) = -a ft -byl, 



(36) 



where a and b are matrices in flavour space, u is the four-velocity of the plasma and k = (u, k) is the 
external momentum. We have neglected the contribution proportional to the masses of the charged 
leptons. In the plasma rest frame and the mass basis, 



Re(£(u;, k))7o = —h(u, k) — a(u>, k)a • k, 



(37) 
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where h(u>, k) = a(u, k)u; + b(u, k) and is given by 



with 



fi 

f\V,a 

fz 



9 - [f z H( mi , M z ) + f H H( mi , M H )}6, t + fw,aH(M a , M w ) 



Qfg 2 s 2 w (L + R) 



' 1 + f 



— (T 3 



2 

A_ 2 

2 



4 



c 



ii' 



i? 



(38) 

(39) 
(40) 

(41) 

(42) 



where L, R are the chiral projectors, Aj = rrii/Mw, mi are the masses of the external flavours, and M Q 
are the masses of the neutral leptons inside the loop. The function H(Mp, Mb) can be found in ||. 
The dispersion relations of the quasi-particles are then given by 



#-Re(E(Jfe)) = 0. 



(43) 



Since these no longer are Lorentz invariant, it is not possible to simplify the calculation of the reflection 
amplitudes by boostings them to the frame where k x ,k y = 0, as we did in section 4. The realistic 
computation in three dimensions then becomes very involved. However, since our main interest is 
to study whether the enhancement found in [18[] is present, and this can already be seen in the one- 
dimensional problem, we restrict our discussion to this simpler case. 

Our objective is to compute the number current of II reflected on the wall, which for small wall 
velocities and using unitarity and CPT, is given by 



j C p = -2 



duj 

T J 2k 1 



u v 2 a /fM[1 - f F (w)]A CP (uj) 



where v g = || = 1/3 is the group velocity, 



A CP (u;) = £(L%- 



\fi. 



(44) 



(45) 



'■j 



with A r the reflection amplitudes on the wall, and is the unboosted Fermi distribution of quasi- 

particles in the plasma rest frame: 

f F (uj) = l/(e^ T + l). (46) 

To lowest order in the wall velocity, the scattering problem can be approximately solved in the 
rest frame of the wall, neglecting the corrections to the dispersion relations of the quasi-particles due 
to the small boost (which are proportional to v w and are negligible at lowest order). In this frame, 
to leading order in T and neglecting the flavour-non-diagonal corrections in fl38|) , the quasi-particles 
propagate according to the effective Hamiltonian 



H eff 



-la z d z + u u R 

\a z d z + uj u l 



+ 0(z) 



-\o z d z + uj h R 



\a z d z + w\ 



(47) 
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where uj^r, satisfying 

^!2 + <M,0)=0, (48) 

are the thermal masses in the unbroken and broken phases respectively (the functions h contain only 
the leading T flavour-diagonal corrections in (0)). This effective Hamiltonian is only valid for low 
momentum compared to the thermal masses u>^}r- Since the reflection of quasi-particles on the wall 

will occur for k z < m << this approximation is justified. 

In order to obtain a non-vanishing CP asymmetry, both the subleading corrections in T (which 
introduce the dependence on the NHL masses) and the flavour-non-diagonal terms (which contain 
the mixings) in (j38|) are needed. After including these corrections, we get the following effective 
Hamiltonian: 



H, 




6(-z)5h R (u°,0)+9(z)5h b R (uj ,0) 

9(-z)5h u L (uj ,0) + 9(z)5h b L (uj ,0) 



(49) 



where Sh = h — h, contain the subleading effects in T and the flavour-non-diagonal electroweak 
corrections. The reflection amplitudes of quasi-particles on the wall can then be obtained by first 
solving for eigenstates of the unperturbed Hamiltonian (fT7|), which are superpositions of incoming, 



reflected and transmitted plane waves, and then perturbing in the extra terms of fl49). 



Up to now, we have neglected the imaginary part of the one-loop self-energy ( pq) . This contribution 
is proportional to the damping rate of the quasi-particles, i.e. their inverse lifetime. There is no 
calculation of the damping rate 7 of leptons in the SM, but from the result for pure SU{2) at zero 
IE7I 



momentum 



we can estimate 7 ~ ayyT, i.e. 7 ~ 1 GeV at T = 100 GeV. In refs. 



it 



was shown that the damping effects for quarks in the SM lead to a sizeable suppression of the CP 
asymmetry, because the lifetime ~ 1/(27) °f the quasi-quarks in the plasma, was much smaller than 
their reflection time on the wall ~ 1/m (for the down quarks, which gave the leading contribution). 
In the present case, we expect that the main effect will come from the reflection of the r lepton and, 
according to the previous rough estimate, the lifetime of the quasi-tau would be of the same order of 
magnitude as its reflection time. In this situation, it is not clear whether the damping will have an 
important effect or not. We will first compute the asymmetry neglecting the damping completely and 
at the end of this section we will estimate its effect. As we will see, it leads to a suppression that varies 
rapidly with the exact value of 7 around the region 7 ~ m T . 

Similarly to the SM case, we find that the first effect in the asymmetry appears at 0{aiyy)- Defining 
r = Aq (the unperturbed reflection amplitude), we get 



I -J 



x Im 



33 



\d--\ 2 



+ 



mj [{r. 



°) 2 -(r" 7 ) 



d/*} & dj* 



■13 



dji 



(50) 



where Sh b are flavour-dependent and di 



ul + ui, — 2uo + 



+ 



" l 3 L ' R "~ ' 2 1 2 

just as in the SM, the effect comes from the interference of the 0(aw) 
there is no effect at leading order in T, because at this order 8h b R = 0. 
Substituting the expressions for 5h b L R , A^p(cj) can be written as 



From eq. (p 
terms in 8h' 



.(2) 
*CP 



U3 



a 



W S lm ( K HiaK* Hja K mii 

i,j a<b 



K Hjb )f{m limj )F{M a ,M b ) 



we see that, 
and Sh R , and 



(51) 
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where 



and 



/K, mj) = XiXjlm <; r 



F(M a , M b 



03 



i | ^ 



+ 



m A( r u) 2 - (r 



If] 



2,dadijdji 



(52) 



(2 + Xl)I(M a ) - 2/(0) 
(2 + X 2 b )I(M b ) - 2/(0)" 



[I(M b ) - I(0)U 
[(I(M a ) - J(0)] WH . 



(53) 



/(M a ) = |i/(M a , Mjy) and the subscript indicates at which value of the function iZ is evaluated. 
Equation ([51]) shows explicitly the GIM cancellation for both external and internal lepton masses. 

Of all the terms in eq. (|5l"D , the one corresponding to the pair of external flavours (fi, r) gives 
the largest contribution, because f{m^m T ) is a few orders of magnitude larger than for the other 
combinations, while the experimental bounds on the mixings are of the same order. We restrict to this 
leading term for which the 'joint' bound described in section 3 is 



\lm(K H2a K* H3a K* H2b K mb ) 



< 5 x 10" 



(54) 



independently of the flavour of the heavy leptons a, b. Thus, if we assume that all the mixings (54) 
are of the same order of magnitude, the size of the various terms in the sum over the heavy flavours 
depends only on the function F(M a , M b ). We consider just one of these terms as a prototype, i.e. 



A 2 J« b = a 2 w Im(K H2a K* H3a K* H2b K mb )f(m 2 ,m 3 )F(M a ,M b 



(55) 



where there is no sum over a, b. Since we will allow the two heavy masses M a , M b to vary arbitrarily, 
the largest value obtained for the integrated asymmetry, considering only (^), is also an upper bound 
for the other terms. Thus the final result will be at most three times larger, if the terms add coherently. 

In fig. 3 we show the contribution to the CP asymmetry, A^p 6 (u;). We have taken the following 
values for the masses at the phase transition temperature (T ~ 100 GeV): m M = 69 MeV, m T = 1.176 
GeV, M w = 50 GeV, and the weak coupling is aw = 0.035. We have fixed the mass of one NHL to 
M a = 80 GeV, and we plot the result for different values of the other NHL mass. 

The peaks are situated in regions where the r lepton reflects completely, while the [i does not. The 
amplitude of the peaks is larger than one would expect from naive power counting, implying that the 
suppression in the charged lepton masses is not at work, as found in \TE\. However, in contrast to what 
was obtained in the SM, the two peaks tend to cancel each other, and there is a big suppression in the 
integrated result, since the Fermi factors in eq. (f44l) are approximately constant. For M b = 140 GeV 
the contribution of this term to the integrated result in (|44]) is ~ 10~ 12 t>„, (which turns out to be of 
the same order as the naive estimate). Whether the peaks come with equal or opposite signs seems to 
be very dependent on the relative position of the thermal masses of the different flavours. In this case 
the thermal masses are almost flavour-independent, while in the down sector of the SM there is a big 
shift in the third family thermal masses compared to the other two, due to the top Yukawa. This is 
why there is no such cancellation in that case. The conclusion is that the enhancement found in |18] 
is rather model-dependent and it seems to require large flavour-dependent thermal corrections. 

Finally, we want to estimate the effect of the damping rate which, as discussed before, is not 
negligible compared to the reflection time of the r lepton. As shown in refs. |5], §], the decoherence 
effects induced by the damping rate can be taken into account by including the imaginary part of the 
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L_L 

16.4 




16.6 16.8 
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17.6 



17.8 
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Figure 3: CP asymmetry, A^p b (u), for zero damping rate and NHL masses M a = 80 GeV and Mb 
110 GeV (dashed), 140 GeV (solid) and 200 GeV (dotted). 
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Figure 4: A^ b (u;) for M a = 80 GeV, M b = 140 GeV and damping rate 7 = .5 GeV (solid), I. GeV 
(dashed) and 1.5 GeV (dotted). 
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self-energy in the effective Hamiltonian p7| ) and solving it for spatially damped waves. Since the exact 
value of the damping rate is not known, we have computed A"cp b (u) for different values of 7, namely 
7 = 0.5, 1., and 1.5 GeV. The result is shown in fig. 4. As is clear from the curves, the suppression 
due to the damping increases rapidly when 7 ~ m T . Without a precise determination of this quantity, 
it is thus impossible to estimate the actual suppression, although it is clear that neglecting this effect 
is not justified. 

To summarize, we have found that the CP asymmetry in the reflection of the charged leptons is at 
most of O(10~ 12 ), even neglecting decoherence effects due to interactions in the plasma. 



6 Baryon Asymmetry 

In this section, we calculate the baryon asymmetry induced by the CP asymmetries computed in the 
previous sections. This is a very difficult problem since a microscopic treatment is no longer possible 
and we have to match somehow the microscopic result of jcp with the thermodynamic treatment of 
transport of the chiral lepton number generated at the wall. Strictly speaking, the two problems, 
reflection and transport, are completely coupled and should be solved at the same time. This how- 
ever implies treating a many-body non-equilibrium quantum system, and some approximations are 
necessary. 

We will consider here only the effect obtained from the reflection of the neutral leptons, since the 
CP asymmetry of the charged leptons computed in the last section is far too small. In the neutral 
sector, we have completely neglected the thermal incoherence effects in the reflection. This has been 
shown to be a very bad approximation when the damping rate is comparable to both the height and/or 
width of the wall fl H |14| . However, this is not the case here as we discussed in section 3. We believe 
that, because of this, reflection can be treated independently of transport []. 

The picture is then that near the wall in the unbroken phase, a local density of v B lepton number 
is generated due to reflection. This local density generates a diffusion current in the plasma and 
decays due to sphaleron processes that take place in the unbroken phase as we go away from the wall, 
generating a baryon number density. This picture is only consistent if the reflected particles have 
enough time to diffuse before the wall catches up. This will be true for small velocities of the wall. In 
this case also the incoming flux of particles in the calculation of jcp can be taken to be the thermal 
one, as we assumed in the previous sections. 

The diffusion equations in the wall rest frame read 

d t n B \ ( D B d 2 z - v w d z - 3T9(-z) -T6(-z) \( n B \ 

d t n L ) " \ -ZT6{-z) D L d 2 z - v w d z - T6(-z) )\n L ) 

where V = 9^^- and Tws = K,(awT) A is the weak sphaleron rate with k a coefficient of 0(1) We 
have made the further approximation that the sphaleron rate in the broken phase is zero. This drastic 
approximation can only be justified if the phase transition is strongly first order; v w is the velocity of 
the wall. The constants Di B are the diffusion coefficients for leptons and quarks respectively. Since 
quarks suffer strong interactions, it is clear that 

D B < D L . (57) 

3 If the damping rate is not small compared to other scales in the problem, we do not think one can separate the 
problems of reflection and transport, and a detailed calculation is much more complicated. 

4 Howcvcr, there is a recent claim that damping effects in the plasma suppress the sphaleron rate to 0{a\ v T A ) J29f. 
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We take the values of the diffusion constants estimated in the SM in ref. ||10|| , namely Dl ~ 110/T and 
D B ~ 6/T. These estimates are obtained from the elastic scattering, t-channel vector boson exchange 
diagrams, which are expected to dominate the scattering process. Yukawa interactions are neglected. 
We have not included here any other possible vl number decay process than the sphalerons. Other 
decays through Higgs interactions are obviously possible, but we find that their rate is smaller than 
the sphaleron rate in the unbroken phase, so we can safely neglect them. 

We look for stationary solutions, i.e. d t n BL = 0. In order to solve the equations for n LB (|56|), 
we need to impose boundary conditions on the densities and their derivatives (diffusion currents). We 
will require that riL tB (— oo) = 0, since there is no asymmetry in the incoming thermalized flux seen by 
the wall. At z — > oo, we require that the solutions be constants. These would be precisely the values 
of L and B in the broken phase that will survive the phase transition. At the interphase z = 0, we 
impose continuity of the diffusion current, 

D L ,Bd z n L)B -v w n L)B \t = (58) 

and the existence of the reflected flux is taken into account in a constraint on the lepton density in 
the unbroken phase near the wall, 

n L \ z=0 - = n° (59) 

where n° = jcp/ (vi) and (vi) is the average velocity of the particles in the reflected flux that we define 

as, 



(Vi) 



I Jijk{ ~2Fr(M u Mj, M k ) + (^S^)^(M,, M h M k ) }§/*.(?)[! - f F (p)} 



It is straightforward to obtain the most general solutions jn| of ( |56D in the approximation, 

3D B rT 



V 2 
w 



< 1, (60) 



which is expected to be of O(10 ). The solution is 



with 



n B = C iau e kiz + C 2 a 21 e k2Z 
n L = C iai2 e^ z + C 2 a 22 e k2Z z<0 
n B = B 

n L = L z > 



, v w 3TD B 



D B \ v. 



2 




(61) 



and 




1 

3TD F 



vl(D L /D B -l) 




vl(D B /D L -l) \ _ 
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Now, the constants C\p and B,L can be determined from (|58| ) and (|59|). In the limit Db "C -Dl the 
result for S is 



5 = (63) 

Although the dependence on the wall velocity seems to have a singular limit when v w —>■ (there is 
only one power of the wall velocity in n°), this is only because we have made the approximation that 
the ratio — ^ < 1, In the limit v w — > 0, this approximation is obviously not valid and indeed the 
solution of the diffusion equations in this case gives B = 0, as expected. 

Finally, in order to compare this result with the experimental one B/s ~ (4-6) x 10~ n , we need 
to divide by the entropy at the temperature of the phase transition, s = ^-g*T 3 , where g* ~ 110 
counts the degrees of freedom of the relativistic particles at the electroweak phase transition. Putting 
everything together we find generically a effect of the order of 

YD 

B/s Ji 23 x 1(T 4 . (64) 



If we assume v(T)/u(T) = v(0)/u(0), we can use the experimental bound J123 < 1CT 7 flltt ). Considering 
the values quoted in the literature for the ratio YDl/v w ~ 10~ 2 (within the SM) |10| |24|] P5 |, we get a 
baryon to entropy ratio two orders of magnitude smaller than required. However, the bounds on J123 
only hold if the ratio of the scalar vevs does not vary with the temperature, which is not necessarily 
true. For instance, a variation by a factor of 2 in the right direction (i.e. a larger ratio at T), increases 
the result by two orders of magnitude. This is because the CP asymmetry goes like 0(D/S) 6 which, 
up to Yukawa couplings, is ~ 0(v(T)/u(T)) 6 . An enhancement due to this effect has been suggested 
in the context of two-Higgs models in []10[ . In order to establish whether this enhancement could take 
place, a detailed study of the scalar potential is required, which is beyond the scope of this paper. 



7 Conclusions 

We have considered the possibility that baryogenesis occurs during the weak phase transition in a 
minimal extension of the Standard Model, which contains extra neutral leptons and conserves total 
lepton number. The leading CP asymmetries come from the reflection of both neutral and charged 
leptons on the bubble wall. Due to the large mean free path of the leptons as compared to the 
typical values of the wall thickness, the calculation is done in the thin wall regime. The CP-violating 
phases come from two sources. For the NHL there are unremovable CP phases due to the non-trivial 
space dependence of the mass matrix inside the bubble wall. The effect turns out to be tree level 
and in agreement with naive estimates. It is only suppressed by the mixing angles. For the charged 
leptons there is no tree-level contribution and the CP- violating phases appear in the one-loop thermal 
corrections to the lepton propagation in the plasma. The naive estimate gives a suppression in the 
charged lepton masses and in loop factors (l/47r), besides that in the mixing angles. The result of 



18| suggests that the suppression on the charged lepton masses could be absent; however, we find 



agreement with the naive estimate. We argued that the effect found in [18| requires a large flavour 



dependence of the leading T thermal corrections, which is not the case in this type of models. 

Using the present constraints on the mixing angles, we obtain that the leading effect comes from the 
neutral sector and gives B/s < 10~ n . Assuming SM estimates for the lepton diffusion constant 
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Dl, the wall velocity v w and the sphaleron rate T, we get B/s < 10~ 13 , which even though the errors 
involved are very large, seems too small to account for the observed baryon asymmetry. However, the 
constraints on the mixing angles apply only if the ratio of the scalar vevs at the temperature of the 
phase transition is the same as today, which is not necessarily true. It would interesting to study a 
realistic scalar potential to determine whether this possibility is realized. 

Finally, we want to comment on other scenarios where the baryon asymmetry is also generated at 
the electroweak phase transition, through lepton reflection on the bubble wall. In ref. |13| the singlet 
majoron model was considered. The CP asymmetry in that case was also due to the reflection of 
neutrinos. However, the relevant phase space was around the mass of the r-neutrino ~ O(10 MeV). 
Although the asymmetry obtained was roughly of the correct order of magnitude, we think that thermal 



corrections to the dispersion relation of the v T in the plasma, which were neglected in |13], should be 



taken into account. In particular, from the calculation of the damping rate of neutrinos in this model 



30fl) it is clear that the typical reflection time of the light neutrinos is much larger than the lifetime 
of the quasi-particles in the plasma. In this situation, we expect a considerable suppression in the CP 
asymmetry. In refs. |10l|i~2||, the reflection of r leptons was considered as the leading contribution to 
the baryon asymmetry in the two-Higgs model, in the thin wall regime. The effects of the damping 
rate have also been neglected in this case. The results for the charged lepton contribution to the 
asymmetry in the present work show that this effect could be important. 
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Appendix A. 

In the calculation of the CP asymmetry in the neutral sector (section 4) we have used the exact 
propagator in the presence of a wall in position space. We give here the expression for (eq. (|30|)) 
in the boosted frame, p x = p y = 0: 

C(0) _ o(0) , rr(O) , q(0) (A i \ 

° — °left + ° across + ° right L ) 

with 

+ Q( Zl - z 2 )e- iE ^~ z 0i^«!L + e(-z 2 )e-^+^ 1 ~ az F mT ° / } (A. 2) 

2 2 h + p J 
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-0(, 1 )e(-, 2 )e-'-e-^l^ (l + (A. 3) 



+ e(,),—.i - + ai + ™ 7 oW_ ( , 4) 



and p' = \J E 2 — m 2 . 

The propagator in position space for the massless left-handed neutrino is 

S^(z 2 , Zl ) = &(z 2 - Zl y E ^\ (A. 5) 

Appendix B. 

I\ a = f dzV(z)e^> = 5h{e \ p ^ (E+P ' )SH ~ 1) (B- 1) 

i& = jf d^i jf ^ 2 \/(^i)e iB2l e^ 22 {G( Zl - z 2 )e ip ^~ Z2) + 9(z 2 - Zl )e ip ^ Z2 - Zl) - e ip ^ Zl+Z2) ) 

L J 

Ig = f°° d Zl r dz 2 dz 3 V(z 1 )V(z 2 )V(z 3 )e lE(zi - Z2+zs) e l ^ Z3 (B. 3) 

JO Jo Jz 2 

{(E + Pj )e( Zl - z 2 )e ip ^- Z2) + (E — pj)[Q(z 2 - Zl )e tp ^ Z2 - Zl) - e ^ {zi+Z2) ]} 

poo roo roo roo poo 

Ig k = / dz x / dz 2 \ dz 3 / dz 4 / cfeg^i) • • • V(z 5 )e iE{zi ' Z2+Z3 - Z4+Z5) e mZ5 

Jo Jo Jz 2 Jo Jz4 

{(E + Pj )e( Zl - z 2 y p ^- Z2) + (E — pj)[Q(z 2 - Zl )e ip ^ Z2 - Zl) - e ^^+^)]| 

{{E + p k )Q(z 3 - z A )e iPkizA ~ Zi) + {E- Pk )[Q{z 4 - z 3 )e ipk{z4 - Z3) - e ipfc(23+24) ]} (B. 4) 



r u = 

i% = 

jijk _ 

1 5b — 



liai-Pi) ~ AM 
Il 3a(-Pi) - tlKPi) 
A 5a 



-Pi) ~ III (Pi) 



(B. 5) 

(B. 6) 

(B. 7) 
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